{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "对airmiles数据进行自相关分析"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEICAYAAABS0fM3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3df5xddX3n8dd7fmcySSa/yS8ShFQBWYKdita2SxEV3V2hXavgWsHFpu5KtfaHglrKUrW03S62XbfbrKL4owrSqmmLRUQp/aGUAaNCWEwEQn5nmGRmksy9c399+sc9E26GmclM7p25c+99Px+P+7jnfM/33PM5d5Lzuef7Ped7FBGYmVnjaqp2AGZmVl1OBGZmDc6JwMyswTkRmJk1OCcCM7MG50RgZtbgnAis4iT9F0nfqODnfVDSJ5PpDZJCUkulPr9RSLpW0j+Vsf7XJV1TyZhsbpDvI7BaImkD8DTQGhG56kZTJCmAjRGxs9qxTEbStcA7I+JnplD3ZuCciHjbTMdl1eczAptV/iVfGeN9j/5u7XQ5EdhpkXSDpB9LOippu6RfKFl2UhNE0pTzbkk7gB2SLpG0R9L7JR2StF/SlZLeIOlHkg5L+mDJ+jdL+vwEcSyS9KnkM/ZK+oik5mTZOZL+QdKgpOck3TnJ/nxZ0oGk7oOSzi9Z9oCkd463f5IeTIq/L+mYpLck5b8iaWeyL1slrS5Z/3xJ9yXLDo7uq6R2SR+XtC95fVxSe7Js9Dv7gKQDwKeT7+VuSZ+XNARcO9n3Mc4+/4mk3ZKGJD0i6WeT8suBDwJvSfbp+2O/B0lNkj4saVfyN/yspEXJstHmu2skPZt89x+a6Lu36nMisNP1Y+BngUXA/wA+L2nVJPWvBC4GzkvmzwA6gDXATcD/A94G/GTyub8j6awpxPEZIAecA1wEvBYYPWj/HvANYDGwFvizST7n68BGYAXwKPCFKWybiPi5ZPLCiOiKiDslXQr8PvBmYBWwC/gSgKQFwDeBvwdWJ3Hfn3zGh4BXAJuAC4GXAx8u2dwZwBJgPbA5KbsCuBvoTmKe7PsY6+FkW0uAvwS+LKkjIv4e+BhwZ7JPF46z7rXJ6+eBFwFdwP8eU+dngBcDrwZuknTuBHFYtUWEX36V/QK2AVck09cC/1SyLIBLS+YvAVJAczK/IKlzcUmdR4Ark+mbgc8n0xuSui3ASmAEmFey3tXAt5PpzwJbgLXT3JfuZBuLkvkHKLatM8n+nVMy/yngD0vmu4BsEvvVwPcm2O6PgTeUzL8OeKbkO8sAHSXLbwYeLJk/1fdxUtzjbP8IxYR20ndesvzE90Axef33kmUvTvaxpeRvtLZk+b8CV1X736lf4798RmCnRdLbJW2TNCBpAHgpsGySVXaPme+PiHwynUreD5YsT1E8gE5mPdAK7C+J4y8o/qoHeD8g4F8lPS7pv06wL82Sbk2auoaAZ5JFk+3PZFZTPAsAICKOAf0Uz37WUTzgn3K9ZHp1yXxfRKTHrFP6vZ7q+ziJpN+S9ETSHDZA8exuqvs8XqyjyXnUgZLpYU7997QqceeSTZuk9RSbcl4NfCci8pK2UTzoTmQmLk/bTfEX8LIY5wqiiDgA/AqApJ8BvinpwXjh1T1vpdjEchnFJLCI4q/j0f05DnSW1D/jFHHto3hQJtn2fGApsDeJ+apTrPd4Mn9mUnZil8ZZp7Rs0u+jVNIf8H6Kf8PHI6IgqXSfT/X3Omkfk1hzFJP52lOsa3OMzwjsdMyneKDoA5D0DopnBLMqIvZT7AP4Y0kLkw7MsyX9+ySuX5I0elA6ksRcGOejFlA8gPZTPOB/bMzybcAvSuqUdA5w3ZjlBym2k4/6IvAOSZuSzt6PAQ9FxDPA3wKrJP160jm8QNLFJet9WNJyScso9p2M20l+Ot/HOPuco/g3bJF0E7BwzD5tkDTRMeKLwPsknSWpi+f7FObEJb02PU4ENm0RsR34Y+A7FA8YFwD/XKVw3g60AdspHuzvpthBC/BTwEOSjgFbgfdGxFPjfMZnKTZt7E0+57tjlt9GsX3+IHAHL+xIvhm4I2mOeXNEfBP4HeCvgP3A2SRnARFxFHgN8J8oNp3soNjhCvARoBf4AfBDip3WH5n6VwFM/n2Uupdih/WPkn1Pc3Iz05eT935Jj46z/u3A54AHKd7XkQZ+bZqx2hzhG8rMzBqczwjMzBqcE4GZWYNzIjAza3BOBGZmDa4m7yNYtmxZbNiwodphmJnVlEceeeS5iFg+trwmE8GGDRvo7e2tdhhmZjVF0q7xyt00ZGbW4JwIzMwanBOBmVmDcyIwM2twTgRmZg2uIolA0u3J4+oem2C5JP1p8ui+H0h6WcmyayTtSF7XVCKe8eQLwf1PHORP79/B/U8cJF/wGEtmZlC5y0c/Q/ExdZ+dYPnrKT4GcCPFxxX+OXCxpCXA7wI9FIcIfkTS1og4UqG4gGIS+OVPPcS23QOkMnnmtTWzaV03n7vuYpqbJhtC38ys/lXkjCAiHgQOT1LlCuCzUfRdoDt5vu3rgPsi4nBy8L8PuLwSMZV64MlDbNs9wHAmTwDDmTzbdg/wwJOHKr0pM7OaM1t9BGs4eazzPUnZROUvIGmzpF5JvX19fdPa+OP7hkhl8ieVpTJ5tu8bmtbnmJnVo5rpLI6ILRHRExE9y5e/4A7pSZ2/eiHz2ppPKpvX1sx5qxdOsIaZWeOYrUSwl+JDu0etTcomKq+oS168gk3rulE+A1GgM+kjuOTF4z7T28ysocxWItgKvD25eugVwGDyfNV7gddKWixpMfDapKyimpvE5667mOU7/obuPf/Mn119kTuKzcwSFblqSNIXgUuAZZL2ULwSqBUgIv4vcA/wBmAnMAy8I1l2WNLvAQ8nH3VLREzW6XzamptE58BTdA48xavPXTkTmzAzq0kVSQQRcfUplgfw7gmW3U7xQdhmZlYFNdNZbGZmM8OJwMyswdXkg2nMbGYUCkEhgkJAEETJSCyj00GMmYdIZmJMXWL8+iTrjK3/fMkcMEkok0UZk643/sLJ13le97xW5rdX/rDtRGA2h0UE2XyQP3GATg7SyXu+ECemR5dHMl1chxcsLxQY97NG17W5q2X5fCeCWpMvBA88eYjH9w1x/uqFXPLiFb5ktcHlC0E2XyBXCLK5AtlCgVw+yOXjxHQ2XzhRJ5f3kdlmnhPBDPFAd/UvIk4crDP5ArnRA3w+ObgXCmRyxfdsPsjlC3jQW5uLnAhmSOlAd3DyQHe+j6F2jOTypDMFUtk8qWyedDbPSO75g76bUqweOBHMkMkGunMimFsKhSCdy5PKPH+wTyUHfz+3whqBE8EMGR3obrgkGXigu+rK5AolB/qTf+H7l701MieCGTI60N13frSfaGqhs73VA93NgoggnX2+KSeVSQ782bw7Xs0m4EQwQ0YHunvlL15HZv4K/vjD7/NVQxU2nMlxLJ076aDvX/dm0+dEMIM80F1lZfMFBlNZBoazDKayZHKFaodkVhecCGzOigiG0jkGkwP/sZFctUMyq0tOBDanpDL54q/+VIahVM5X7ZjNAicCq6rcaHNPqvirfyTr5h6z2eZEYLMqIjg6cnJzjzt3zarLicBmXDqbP9HJO5TO+jJOszmmUo+qvBz4E6AZ+GRE3Dpm+W3AzyezncCKiOhOluWBHybLno2IN1YiJquefCEYTJp6BoYzpN3cYzanlZ0IJDUDnwBeA+wBHpa0NSK2j9aJiPeV1P814KKSj0hFxKZy47Ciao14ms7mee7YCAPDbu4xqzWVOCN4ObAzIp4CkPQl4Apg+wT1r6b4cHursGqMeHo0nWX/YJrDxzM++JvVqEo8qnINsLtkfk9S9gKS1gNnAd8qKe6Q1Cvpu5KunGgjkjYn9Xr7+voqEHb9KR3xNDh5xNNKigj6j43w2N5BHts7RP8xJwGzWjbbzyy+Crg7IkqH5VwfET3AW4GPSzp7vBUjYktE9EREz/Lly2cj1poz2YinlZAvBPsHU3xv9wA/OniMo2nf4GVWDyrRNLQXWFcyvzYpG89VwLtLCyJib/L+lKQHKPYf/LgCcTWcmRrxdCSX58BgmkNHR3zFj1kdqsQZwcPARklnSWqjeLDfOraSpJcAi4HvlJQtltSeTC8DXsXEfQt2CqMjniqfgSjQmfQRnO6Ip8dGcuw8dJTvPTvAvoG0k4BZnSr7jCAicpKuB+6lePno7RHxuKRbgN6IGE0KVwFfijipNflc4C8kFSgmpVtLrzay6anUiKeHj2fYP5hiKOWmH7NGUJH7CCLiHuCeMWU3jZm/eZz1/gW4oBIxWNHpjnhaKAR9x0bYP5h+QT+DmdU331nc4DK5AgeH0hwcSpN1049ZQ3IiaFDDmRz7BtL0HxvBA3yaNTYnggYzMJxh/2CageFstUMxsznCiaABFArBc0n7/7Db/81sDCeCOpbNP9/+n8nNXPtPoRBs2z3AM/3H2bB0PpvWddPkZzOb1QwngjpUiCCTK/DoriMz3v5fKAQf+/oT7Dx0jEyuQFtLE+es6OKDrz/XycCsRsz2EBM2g/KF4Nn+YY6N5MjkCrPSCbxt9wA7Dx1jJFcggJFcgZ2HjrFt98DMb9zMKsKJoA5EBAeH0mzbfYS9AymYxauAnuk/TiZ38vMGMrkCz/Qfn70gzKwsbhqqcQPDGXb1D1etE3jD0vm0tTQxUpIM2lqa2LB0flXiMbPp8xlBjRrO5Hhi/xBP7D9a1SuBNq3r5pwVXZArjm/UnvQRbFrXXbWYzGx6nAhqTCZX4Km+Y/xgz+CcuBegqUl88PXn0rX9q8x7+h95z6Ub3VFsVmPcNFQjCoVg32CKfQNp8nPsVuCmJtHWvxP6d/Ky9R+odjhmNk1OBDWg7+gIu48MM+KHwJvZDHAimMOG0lme7R/2k8DMbEY5EcxB6WyeZw8P038sU+1QzKwBOBHMIbl8gb0DKQ4Mpj0iqJnNGieCOaB4Q9gIe44M+5kAZjbrKnL5qKTLJT0paaekG8ZZfq2kPknbktc7S5ZdI2lH8rqmEvHUksPHM3x/zyBPP3fcScDMqqLsMwJJzcAngNcAe4CHJW0d59nDd0bE9WPWXQL8LtBDcWCER5J1j5Qb11x3fCTHrv5hBlPVvxegFnnEU7PKqUTT0MuBnRHxFICkLwFXAFN5CP3rgPsi4nCy7n3A5cAXKxDXnDSSy7P7cIrnjo0QPgE4LR7x1KyyKtE0tAbYXTK/Jykb6z9L+oGkuyWtm+a6SNosqVdSb19fXwXCnl35QrD78DDf3z1I31EngXJ4xFOzypqtISb+BtgQEf8OuA+4Y7ofEBFbIqInInqWL19e8QBnUjZfYNvuAfYcSc25u4JrkUc8NausSiSCvcC6kvm1SdkJEdEfESPJ7CeBn5zqurVsMJXl+EiOVCb/ggOXnb7REU9LecRTs9NXiUTwMLBR0lmS2oCrgK2lFSStKpl9I/BEMn0v8FpJiyUtBl6blNW0VCbPkweOsn3fkM8AZoBHPDWrrLI7iyMiJ+l6igfwZuD2iHhc0i1Ab0RsBd4j6Y1ADjgMXJuse1jS71FMJgC3jHYc16JsvsCeIykODqXdBzCDRkc8/dX3/ib5rpVc/67NvmrIrAwVuaEsIu4B7hlTdlPJ9I3AjROseztweyXiqJZCITgwlGbvQIqc7wWYFR7x1KxyfGdxmfqPjfDs4WHSHhnUzGqUE8FpOjaS45nnjntkUDOreU4E01S8IWyYvqMeGdTM6oMTwRTlC8G+gRT7BlIeGdTM6ooTwRQdGEqz50iq2mGYmVWcH15vZtbgfEZgNkUe8dTqlROB2RR4xFOrZ24aMpsCj3hq9cyJwGwKPOKp1TMnArMp8IinVs+cCMymwCOeWj1zIjCbgtERT7u2f5V5T/8j77l0ozuKrW44EZhN0eiIp/N2/TMvW7/YScDqhhOBmVmDcyIwM2twFUkEki6X9KSknZJuGGf5b0jaLukHku6XtL5kWV7StuS1dey6ZmY2s8q+s1hSM/AJ4DXAHuBhSVsjYntJte8BPRExLOm/AX8IvCVZloqITeXGYWZmp6cSZwQvB3ZGxFMRkQG+BFxRWiEivh0Rw8nsd4G1FdiuWd0rFIJHdx3hrx/dw6O7jlDwGOg2Ayox1tAaYHfJ/B7g4knqXwd8vWS+Q1IvxQfb3xoRXx1vJUmbgc0AZ555ZlkBm9UCj29ks2VWO4slvQ3oAf6opHh9RPQAbwU+Luns8daNiC0R0RMRPcuXL5+FaM2qy+Mb2WypRCLYC6wrmV+blJ1E0mXAh4A3RsTIaHlE7E3enwIeAC6qQExmNc/jG9lsqUQieBjYKOksSW3AVcBJV/9Iugj4C4pJ4FBJ+WJJ7cn0MuBVQGkns1nD8vhGNlvKTgQRkQOuB+4FngDuiojHJd0i6Y1JtT8CuoAvj7lM9FygV9L3gW9T7CNwIjDD4xvZ7KnIg2ki4h7gnjFlN5VMXzbBev8CXFCJGMzqzej4Rr/63t8k37WS69+12U9FsxnhO4vN5jCPb2SzwYnAzKzBORGYmTU4JwIzswbnRGBm1uAqctWQmdWPQiHYtnuAZ/qPs2HpfF+p1ACcCMzsBI9v1JjcNGRmJ3h8o8bkRGBmJ3h8o8bkRGBmJ3h8o8bkRGBmJ3h8o8bkRGBmJ4yOb9S1/avMe/ofec+lG91R3ACcCMzsJB7fqPE4EZiZNTgnAjOzBucbysys6nw3c3U5EZhZVflu5uqrSNOQpMslPSlpp6QbxlneLunOZPlDkjaULLsxKX9S0usqEY+Z1Q7fzVx9ZScCSc3AJ4DXA+cBV0s6b0y164AjEXEOcBvwB8m651F82P35wOXA/0k+z8wahO9mrj5FRHkfIL0SuDkiXpfM3wgQEb9fUufepM53JLUAB4DlwA2ldUvrTbbNJevPjdd88PZpx7rt+9sA2HThpmmvm8kVSOfy015vx/bHANh43kunve7p8ja9zVra5tF0jr0DKUoPRRKs6Z7Hgg63Xpea19pMa/Pp/36/610//UhE9Iwtr8S3vAbYXTK/B7h4ojoRkZM0CCxNyr87Zt01421E0mZgM0DXqrNPK9DTSQDlms3/vN6mt1mL2+xqb2ZeazOpbJ6IYhKY19pMV/vMNw7Ue5KdqppJtxGxBdgC0NPTE3f+6itndft7B1I82z88q9s0axTVumro3X95IwA3vX/rjG+rEtt80fL5rFzYcdrbvutd45dXIhHsBdaVzK9NysarsydpGloE9E9xXTOrc01N4mXrF/Oy9YurHUpDqsRVQw8DGyWdJamNYufv2FS3FbgmmX4T8K0odk5sBa5Krio6C9gI/GsFYjIzsykq+4wgafO/HrgXaAZuj4jHJd0C9EbEVuBTwOck7QQOU0wWJPXuArYDOeDdETH9HlkzMzttFekjiIh7gHvGlN1UMp0GfmmCdT8KfLQScZiZTVWhEGSWnkO+ayWP7jrS0Hcz10xnsZlZpYzezXzsvCuhuYU//daOhr6b2YPOmVnDGb2bmZY2UFPD383sRGBmDcd3M5/MicDMGo6fzXwyJwIzazijz2Zub2lC0PDPZnZnsZk1nNFnM/sZCEVOBGbWkKpxN/NcvWTVicDMbBbM5UtW3UdgZjYL5vIlq04EZmazYC5fsupEYGY2C+byJatOBGZms2AuX7LqzmIzs1kwly9ZdSIwM5slc/UBPG4aMjNrcE4EZmYNzonAzKzBlZUIJC2RdJ+kHcn7Cxq+JG2S9B1Jj0v6gaS3lCz7jKSnJW1LXpvKicfMzKav3DOCG4D7I2IjcH8yP9Yw8PaIOB+4HPi4pNLrpX47IjYlr21lxmNmZtNUbiK4Argjmb4DuHJshYj4UUTsSKb3AYeA5WVu18zMKqTcRLAyIvYn0weAlZNVlvRyoA34cUnxR5Mmo9sktU+y7mZJvZJ6+/r6ygzbzMxGnTIRSPqmpMfGeV1RWi8iAohJPmcV8DngHRExOuDGjcBLgJ8ClgAfmGj9iNgSET0R0bN8uU8ozMwq5ZQ3lEXEZRMtk3RQ0qqI2J8c6A9NUG8h8HfAhyLiuyWfPXo2MSLp08BvTSt6MzMrW7lNQ1uBa5Lpa4Cvja0gqQ34CvDZiLh7zLJVybso9i88VmY8ZmY2TeUmgluB10jaAVyWzCOpR9InkzpvBn4OuHacy0S/IOmHwA+BZcBHyozHzMymqayxhiKiH3j1OOW9wDuT6c8Dn59g/UvL2b6ZmZXPdxabmTU4JwIzswbnRGBm1uCcCMzMGpwTgZlZg3MiMDNrcE4EZmYNzonAzKzBORGYmTU4J4Ip6mproa1F1Q7DzKzinAimaFFnKxeu7WZN9zyanA/MrI44EUxDS3MTZy7t5MJ13Szraqt2OGZmFeFEcBo6WpvZuHIBL12zkAUdZY3bZ2ZWdU4EZVjQ0cpL1yziJ1Z20d7qr9LMapN/zlbA0q52Fne2cWAozd6BFLn8hE/sNDObc/wztkKamsTq7nlsWtfNGYs6kDuUzaxGOBFUWGtzE2ctm8+Fa7tZPL+12uGYmZ1SWYlA0hJJ90nakbwvnqBevuQxlVtLys+S9JCknZLuTJ5vXBfmtTXzkjMWct6qhcxvb652OGZmEyr3jOAG4P6I2Ajcn8yPJxURm5LXG0vK/wC4LSLOAY4A15UZz5yzqLOVC9Ys4uzl831DmpnNSeUmgiuAO5LpO4Arp7qiJAGXAnefzvq1RBIrFnawad1i1i6eR7PvSDOzOaTcRLAyIvYn0weAlRPU65DUK+m7kkYP9kuBgYjIJfN7gDUTbUjS5uQzevv6+soMuzqam8S6JZ1cuG4Ryxe0VzscMzNgCpePSvomcMY4iz5UOhMRIWmi6ybXR8ReSS8CviXph8DgdAKNiC3AFoCenp6avj6zvaWZc1Z0ccaiDnb1H2colTv1SmZmM+SUiSAiLptomaSDklZFxH5Jq4BDE3zG3uT9KUkPABcBfwV0S2pJzgrWAntPYx9qVld7C+evXsTh4xl29R8nnS1UOyQza0DlNg1tBa5Jpq8Bvja2gqTFktqT6WXAq4DtERHAt4E3TbZ+I1gyv41N67rZsKyTlmb3H5jZ7Co3EdwKvEbSDuCyZB5JPZI+mdQ5F+iV9H2KB/5bI2J7suwDwG9I2kmxz+BTZcZTsySxatE8LlrXzSrfkGZms0jFH+a1paenJ3p7e6sdxoxKZ/Ps6h/m8PFMtUMxszniRcvns3Jhx2mvL+mRiOgZW+47i+eojtZmXnzGAs5bvZCudg8JZWYzx4lgjls0r5UL1i7inBVdtLX4z2VmleefmjVi+YJ2ls5vY99gin0DafKF2mvSM7O5yT8xa0hTk1i7uJNN67pZsbDdHcpmVhFOBDWoraWJs5d3ccGaRSya5xFOzaw8TgQ1bH57C+etXshLzljAvDaPcGpmp8d9BHVg8fw2ujtbOTg0wp4jw2T9hDQzmwYngjohiTMWdbCsq419A2n2D6Zwf7KZTYUTQZ1paW7izKWdrOru4OBQmoNDaTI5ZwQzm5gTQZ1qbW5i7eJOVi+ax3PHRtg/mGY4k692WGY2BzkR1LmmpuJDcVYs7GBgOMO+gTSDqWy1wzKzOcSJoIF0d7bR3dnGcCbHvoE0/cdG3I9gZk4EjaizrYVzVnRx5pLOE/0IvtLIrHE5ETSwtpYm1i3pZHV3sR9h30DKD8cxa0BOBEZzk1i5sIMVC9o5Mpxl/2DKj880ayBOBHaCJJbMb2PJ/DaOjeTYP5Ci/3iGGnxkhZlNQ1lDTEhaIuk+STuS98Xj1Pl5SdtKXmlJVybLPiPp6ZJlm8qJxyqnq72FjSsXcNGZ3azu7vAjNM3qWLljDd0A3B8RG4H7k/mTRMS3I2JTRGwCLgWGgW+UVPnt0eURsa3MeKzC2luaWb90Pi87czEblnXS3urhqczqTbn/q68A7kim7wCuPEX9NwFfj4jhMrdrs6y56flnKv/Eyi4WdLhV0axelJsIVkbE/mT6ALDyFPWvAr44puyjkn4g6TZJ7ROtKGmzpF5JvX19fWWEbOWQxNKudl66ZhHnr1nI0q42PxfBrMad8medpG8CZ4yz6EOlMxERkibsVpS0CrgAuLek+EaKCaQN2AJ8ALhlvPUjYktSh56eHndfzgELO1pZ2NFKOptn/2DxBjXfj2BWe06ZCCLisomWSTooaVVE7E8O9Icm+ag3A1+JiBPjG5ScTYxI+jTwW1OM2+aQjtZmzlo2nw1LOzmeyTMwnGEwleVoOucrjsxqQLkNvVuBa4Bbk/evTVL3aopnACeUJBFR7F94rMx4rIok0dXeQld7C2sXQ74QDKayDKayDAxnfLOa2RxVbiK4FbhL0nXALoq/+pHUA7wrIt6ZzG8A1gH/MGb9L0haDgjYBryrzHhsDmluev6+BJhPOps/kRgGU1lybkYymxPKSgQR0Q+8epzyXuCdJfPPAGvGqXdpOdu32tLR2kxHazMrF3YQERwdyTE4XEwKx0bcjGRWLb4G0KpC0onO5nVALl94vhkplWXEzUhms8aJwOaEluYmlna1s7SreAVxKpNPkkKGoVSOvMfLNpsxTgQ2J81ra2ZeWzNnLCo2Iw2lcwylsgwMZzmecTOSWSU5EdicJ4lF81pZNK+VdUsgmzQjDQxnOT6SI53N+wE7ZmVwIrCa09rcxLKudpYlzUgRwUiuQCqTJ5VNXpk86WzeN7iZTYETgdU8SSeuSBo7/G02XyCVzZPO5ElnCycSRTqbd/OSWcKJwOpaa3MTrc1NLOxoPak8Ik5KDKNnEKls3vc3WMNxIrCGJOlEh/RYmVzhxFlDKpMnnSu+j+QKPouwuuREYDZGW0sTbS1NLJp38llEoRCkc3lGsgWyhQK5fJDLx4npbL5ANl8gVwifVVhNcSIwm6KmJtHZ1kJn26nrRgTZfJArFMgmSWI0WRQTRaFkeXHaZxtWLU4EZjNAEm0tom0aj/zIJUniRNJIkkguXyBfCAoBhYjiq1CcjqQsH0FEUiepazZVTgRmc0RLcxMtzcUxmSqhmBBGk8RoEkkSSOH5BFIYu7zw/DxAUDzDGZ0GSs5enj+TGbssSpeVlJ28/guXzbbJzsQmiygmWXHy9YrtWIsAAASVSURBVCaNZtL1ZuohUE4EZnWqqUk0If8nt1Pyk8jNzBqcE4GZWYNzIjAza3BOBGZmDa6sRCDplyQ9LqmQPJ5yonqXS3pS0k5JN5SUnyXpoaT8TklTuELbzMwqqdwzgseAXwQenKiCpGbgE8DrgfOAqyWdlyz+A+C2iDgHOAJcV2Y8ZmY2TWUlgoh4IiKePEW1lwM7I+KpiMgAXwKukCTgUuDupN4dwJXlxGNmZtM3G30Ea4DdJfN7krKlwEBE5MaUj0vSZkm9knr7+vpmLFgzs0ZzyntNJH0TOGOcRR+KiK9VPqTxRcQWYAtAT0+Pb6A3M6uQUyaCiLiszG3sBdaVzK9NyvqBbkktyVnBaPkpPfLII89J2nWa8SwDnjvNdWuJ97O+eD/rS7X2c/14hbNx9/nDwEZJZ1E80F8FvDUiQtK3gTdR7De4BpjSGUZELD/dYCT1RsSEVzjVC+9nffF+1pe5tp/lXj76C5L2AK8E/k7SvUn5akn3ACS/9q8H7gWeAO6KiMeTj/gA8BuSdlLsM/hUOfGYmdn0lXVGEBFfAb4yTvk+4A0l8/cA94xT7ymKVxWZmVmVNOKdxVuqHcAs8X7WF+9nfZlT+6nJxtQ2M7P614hnBGZmVsKJwMyswTVUIpho8Lt6ImmdpG9L2p4MCPjeasc0UyQ1S/qepL+tdiwzSVK3pLsl/X9JT0h6ZbVjqjRJ70v+vT4m6YuSOqodU6VIul3SIUmPlZQtkXSfpB3J++JqxtgwieAUg9/VkxzwmxFxHvAK4N11up8A76V4SXK9+xPg7yPiJcCF1Nk+S1oDvAfoiYiXAs0U7zeqF58BLh9TdgNwf0RsBO5P5qumYRIBEwx+V+WYKi4i9kfEo8n0UYoHjQnHcKpVktYC/wH4ZLVjmUmSFgE/R3KPTURkImKgulHNiBZgnqQWoBPYV+V4KiYiHgQOjym+guJAmzAHBtxspEQw0eB3dUvSBuAi4KHqRjIjPg68HyhUO5AZdhbQB3w6aQb7pKT51Q6qkiJiL/A/gWeB/cBgRHyjulHNuJURsT+ZPgCsrGYwjZQIGoqkLuCvgF+PiKFqx1NJkv4jcCgiHql2LLOgBXgZ8OcRcRFwnCo3I1Ra0j5+BcWktxqYL+lt1Y1q9kTxGv6qXsffSIlgosHv6o6kVopJ4AsR8dfVjmcGvAp4o6RnKDbxXSrp89UNacbsAfZExOhZ3d0UE0M9uQx4OiL6IiIL/DXw01WOaaYdlLQKIHk/VM1gGikRnBj8Lnkk5lXA1irHVHHJA38+BTwREf+r2vHMhIi4MSLWRsQGin/Hb0VEXf6CjIgDwG5JL06KXg1sr2JIM+FZ4BWSOpN/v6+mzjrEx7GV4kCbMI0BN2fKbIw+OidERE7S6OB3zcDtJYPf1ZNXAb8M/FDStqTsg8l4T1abfg34QvID5ingHVWOp6Ii4iFJdwOPUrzq7XvMsSEYyiHpi8AlwLJkkM7fBW4F7pJ0HbALeHP1IvQQE2ZmDa+RmobMzGwcTgRmZg3OicDMrME5EZiZNTgnAjOzBudEYGbW4JwIzMwa3L8BpvJKhJkZ07cAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEICAYAAABS0fM3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3df5xddX3n8dd7fmcySSa/yS8ShFQBWYKdita2SxEV3V2hXavgWsHFpu5KtfaHglrKUrW03S62XbfbrKL4owrSqmmLRUQp/aGUAaNCWEwEQn5nmGRmksy9c399+sc9E26GmclM7p25c+99Px+P+7jnfM/33PM5d5Lzuef7Ped7FBGYmVnjaqp2AGZmVl1OBGZmDc6JwMyswTkRmJk1OCcCM7MG50RgZtbgnAis4iT9F0nfqODnfVDSJ5PpDZJCUkulPr9RSLpW0j+Vsf7XJV1TyZhsbpDvI7BaImkD8DTQGhG56kZTJCmAjRGxs9qxTEbStcA7I+JnplD3ZuCciHjbTMdl1eczAptV/iVfGeN9j/5u7XQ5EdhpkXSDpB9LOippu6RfKFl2UhNE0pTzbkk7gB2SLpG0R9L7JR2StF/SlZLeIOlHkg5L+mDJ+jdL+vwEcSyS9KnkM/ZK+oik5mTZOZL+QdKgpOck3TnJ/nxZ0oGk7oOSzi9Z9oCkd463f5IeTIq/L+mYpLck5b8iaWeyL1slrS5Z/3xJ9yXLDo7uq6R2SR+XtC95fVxSe7Js9Dv7gKQDwKeT7+VuSZ+XNARcO9n3Mc4+/4mk3ZKGJD0i6WeT8suBDwJvSfbp+2O/B0lNkj4saVfyN/yspEXJstHmu2skPZt89x+a6Lu36nMisNP1Y+BngUXA/wA+L2nVJPWvBC4GzkvmzwA6gDXATcD/A94G/GTyub8j6awpxPEZIAecA1wEvBYYPWj/HvANYDGwFvizST7n68BGYAXwKPCFKWybiPi5ZPLCiOiKiDslXQr8PvBmYBWwC/gSgKQFwDeBvwdWJ3Hfn3zGh4BXAJuAC4GXAx8u2dwZwBJgPbA5KbsCuBvoTmKe7PsY6+FkW0uAvwS+LKkjIv4e+BhwZ7JPF46z7rXJ6+eBFwFdwP8eU+dngBcDrwZuknTuBHFYtUWEX36V/QK2AVck09cC/1SyLIBLS+YvAVJAczK/IKlzcUmdR4Ark+mbgc8n0xuSui3ASmAEmFey3tXAt5PpzwJbgLXT3JfuZBuLkvkHKLatM8n+nVMy/yngD0vmu4BsEvvVwPcm2O6PgTeUzL8OeKbkO8sAHSXLbwYeLJk/1fdxUtzjbP8IxYR20ndesvzE90Axef33kmUvTvaxpeRvtLZk+b8CV1X736lf4798RmCnRdLbJW2TNCBpAHgpsGySVXaPme+PiHwynUreD5YsT1E8gE5mPdAK7C+J4y8o/qoHeD8g4F8lPS7pv06wL82Sbk2auoaAZ5JFk+3PZFZTPAsAICKOAf0Uz37WUTzgn3K9ZHp1yXxfRKTHrFP6vZ7q+ziJpN+S9ETSHDZA8exuqvs8XqyjyXnUgZLpYU7997QqceeSTZuk9RSbcl4NfCci8pK2UTzoTmQmLk/bTfEX8LIY5wqiiDgA/AqApJ8BvinpwXjh1T1vpdjEchnFJLCI4q/j0f05DnSW1D/jFHHto3hQJtn2fGApsDeJ+apTrPd4Mn9mUnZil8ZZp7Rs0u+jVNIf8H6Kf8PHI6IgqXSfT/X3Omkfk1hzFJP52lOsa3OMzwjsdMyneKDoA5D0DopnBLMqIvZT7AP4Y0kLkw7MsyX9+ySuX5I0elA6ksRcGOejFlA8gPZTPOB/bMzybcAvSuqUdA5w3ZjlBym2k4/6IvAOSZuSzt6PAQ9FxDPA3wKrJP160jm8QNLFJet9WNJyScso9p2M20l+Ot/HOPuco/g3bJF0E7BwzD5tkDTRMeKLwPsknSWpi+f7FObEJb02PU4ENm0RsR34Y+A7FA8YFwD/XKVw3g60AdspHuzvpthBC/BTwEOSjgFbgfdGxFPjfMZnKTZt7E0+57tjlt9GsX3+IHAHL+xIvhm4I2mOeXNEfBP4HeCvgP3A2SRnARFxFHgN8J8oNp3soNjhCvARoBf4AfBDip3WH5n6VwFM/n2Uupdih/WPkn1Pc3Iz05eT935Jj46z/u3A54AHKd7XkQZ+bZqx2hzhG8rMzBqczwjMzBqcE4GZWYNzIjAza3BOBGZmDa4m7yNYtmxZbNiwodphmJnVlEceeeS5iFg+trwmE8GGDRvo7e2tdhhmZjVF0q7xyt00ZGbW4JwIzMwanBOBmVmDcyIwM2twTgRmZg2uIolA0u3J4+oem2C5JP1p8ui+H0h6WcmyayTtSF7XVCKe8eQLwf1PHORP79/B/U8cJF/wGEtmZlC5y0c/Q/ExdZ+dYPnrKT4GcCPFxxX+OXCxpCXA7wI9FIcIfkTS1og4UqG4gGIS+OVPPcS23QOkMnnmtTWzaV03n7vuYpqbJhtC38ys/lXkjCAiHgQOT1LlCuCzUfRdoDt5vu3rgPsi4nBy8L8PuLwSMZV64MlDbNs9wHAmTwDDmTzbdg/wwJOHKr0pM7OaM1t9BGs4eazzPUnZROUvIGmzpF5JvX19fdPa+OP7hkhl8ieVpTJ5tu8bmtbnmJnVo5rpLI6ILRHRExE9y5e/4A7pSZ2/eiHz2ppPKpvX1sx5qxdOsIaZWeOYrUSwl+JDu0etTcomKq+oS168gk3rulE+A1GgM+kjuOTF4z7T28ysocxWItgKvD25eugVwGDyfNV7gddKWixpMfDapKyimpvE5667mOU7/obuPf/Mn119kTuKzcwSFblqSNIXgUuAZZL2ULwSqBUgIv4vcA/wBmAnMAy8I1l2WNLvAQ8nH3VLREzW6XzamptE58BTdA48xavPXTkTmzAzq0kVSQQRcfUplgfw7gmW3U7xQdhmZlYFNdNZbGZmM8OJwMyswdXkg2nMbGYUCkEhgkJAEETJSCyj00GMmYdIZmJMXWL8+iTrjK3/fMkcMEkok0UZk643/sLJ13le97xW5rdX/rDtRGA2h0UE2XyQP3GATg7SyXu+ECemR5dHMl1chxcsLxQY97NG17W5q2X5fCeCWpMvBA88eYjH9w1x/uqFXPLiFb5ktcHlC0E2XyBXCLK5AtlCgVw+yOXjxHQ2XzhRJ5f3kdlmnhPBDPFAd/UvIk4crDP5ArnRA3w+ObgXCmRyxfdsPsjlC3jQW5uLnAhmSOlAd3DyQHe+j6F2jOTypDMFUtk8qWyedDbPSO75g76bUqweOBHMkMkGunMimFsKhSCdy5PKPH+wTyUHfz+3whqBE8EMGR3obrgkGXigu+rK5AolB/qTf+H7l701MieCGTI60N13frSfaGqhs73VA93NgoggnX2+KSeVSQ782bw7Xs0m4EQwQ0YHunvlL15HZv4K/vjD7/NVQxU2nMlxLJ076aDvX/dm0+dEMIM80F1lZfMFBlNZBoazDKayZHKFaodkVhecCGzOigiG0jkGkwP/sZFctUMyq0tOBDanpDL54q/+VIahVM5X7ZjNAicCq6rcaHNPqvirfyTr5h6z2eZEYLMqIjg6cnJzjzt3zarLicBmXDqbP9HJO5TO+jJOszmmUo+qvBz4E6AZ+GRE3Dpm+W3AzyezncCKiOhOluWBHybLno2IN1YiJquefCEYTJp6BoYzpN3cYzanlZ0IJDUDnwBeA+wBHpa0NSK2j9aJiPeV1P814KKSj0hFxKZy47Ciao14ms7mee7YCAPDbu4xqzWVOCN4ObAzIp4CkPQl4Apg+wT1r6b4cHursGqMeHo0nWX/YJrDxzM++JvVqEo8qnINsLtkfk9S9gKS1gNnAd8qKe6Q1Cvpu5KunGgjkjYn9Xr7+voqEHb9KR3xNDh5xNNKigj6j43w2N5BHts7RP8xJwGzWjbbzyy+Crg7IkqH5VwfET3AW4GPSzp7vBUjYktE9EREz/Lly2cj1poz2YinlZAvBPsHU3xv9wA/OniMo2nf4GVWDyrRNLQXWFcyvzYpG89VwLtLCyJib/L+lKQHKPYf/LgCcTWcmRrxdCSX58BgmkNHR3zFj1kdqsQZwcPARklnSWqjeLDfOraSpJcAi4HvlJQtltSeTC8DXsXEfQt2CqMjniqfgSjQmfQRnO6Ip8dGcuw8dJTvPTvAvoG0k4BZnSr7jCAicpKuB+6lePno7RHxuKRbgN6IGE0KVwFfijipNflc4C8kFSgmpVtLrzay6anUiKeHj2fYP5hiKOWmH7NGUJH7CCLiHuCeMWU3jZm/eZz1/gW4oBIxWNHpjnhaKAR9x0bYP5h+QT+DmdU331nc4DK5AgeH0hwcSpN1049ZQ3IiaFDDmRz7BtL0HxvBA3yaNTYnggYzMJxh/2CageFstUMxsznCiaABFArBc0n7/7Db/81sDCeCOpbNP9/+n8nNXPtPoRBs2z3AM/3H2bB0PpvWddPkZzOb1QwngjpUiCCTK/DoriMz3v5fKAQf+/oT7Dx0jEyuQFtLE+es6OKDrz/XycCsRsz2EBM2g/KF4Nn+YY6N5MjkCrPSCbxt9wA7Dx1jJFcggJFcgZ2HjrFt98DMb9zMKsKJoA5EBAeH0mzbfYS9AymYxauAnuk/TiZ38vMGMrkCz/Qfn70gzKwsbhqqcQPDGXb1D1etE3jD0vm0tTQxUpIM2lqa2LB0flXiMbPp8xlBjRrO5Hhi/xBP7D9a1SuBNq3r5pwVXZArjm/UnvQRbFrXXbWYzGx6nAhqTCZX4Km+Y/xgz+CcuBegqUl88PXn0rX9q8x7+h95z6Ub3VFsVmPcNFQjCoVg32CKfQNp8nPsVuCmJtHWvxP6d/Ky9R+odjhmNk1OBDWg7+gIu48MM+KHwJvZDHAimMOG0lme7R/2k8DMbEY5EcxB6WyeZw8P038sU+1QzKwBOBHMIbl8gb0DKQ4Mpj0iqJnNGieCOaB4Q9gIe44M+5kAZjbrKnL5qKTLJT0paaekG8ZZfq2kPknbktc7S5ZdI2lH8rqmEvHUksPHM3x/zyBPP3fcScDMqqLsMwJJzcAngNcAe4CHJW0d59nDd0bE9WPWXQL8LtBDcWCER5J1j5Qb11x3fCTHrv5hBlPVvxegFnnEU7PKqUTT0MuBnRHxFICkLwFXAFN5CP3rgPsi4nCy7n3A5cAXKxDXnDSSy7P7cIrnjo0QPgE4LR7x1KyyKtE0tAbYXTK/Jykb6z9L+oGkuyWtm+a6SNosqVdSb19fXwXCnl35QrD78DDf3z1I31EngXJ4xFOzypqtISb+BtgQEf8OuA+4Y7ofEBFbIqInInqWL19e8QBnUjZfYNvuAfYcSc25u4JrkUc8NausSiSCvcC6kvm1SdkJEdEfESPJ7CeBn5zqurVsMJXl+EiOVCb/ggOXnb7REU9LecRTs9NXiUTwMLBR0lmS2oCrgK2lFSStKpl9I/BEMn0v8FpJiyUtBl6blNW0VCbPkweOsn3fkM8AZoBHPDWrrLI7iyMiJ+l6igfwZuD2iHhc0i1Ab0RsBd4j6Y1ADjgMXJuse1jS71FMJgC3jHYc16JsvsCeIykODqXdBzCDRkc8/dX3/ib5rpVc/67NvmrIrAwVuaEsIu4B7hlTdlPJ9I3AjROseztweyXiqJZCITgwlGbvQIqc7wWYFR7x1KxyfGdxmfqPjfDs4WHSHhnUzGqUE8FpOjaS45nnjntkUDOreU4E01S8IWyYvqMeGdTM6oMTwRTlC8G+gRT7BlIeGdTM6ooTwRQdGEqz50iq2mGYmVWcH15vZtbgfEZgNkUe8dTqlROB2RR4xFOrZ24aMpsCj3hq9cyJwGwKPOKp1TMnArMp8IinVs+cCMymwCOeWj1zIjCbgtERT7u2f5V5T/8j77l0ozuKrW44EZhN0eiIp/N2/TMvW7/YScDqhhOBmVmDcyIwM2twFUkEki6X9KSknZJuGGf5b0jaLukHku6XtL5kWV7StuS1dey6ZmY2s8q+s1hSM/AJ4DXAHuBhSVsjYntJte8BPRExLOm/AX8IvCVZloqITeXGYWZmp6cSZwQvB3ZGxFMRkQG+BFxRWiEivh0Rw8nsd4G1FdiuWd0rFIJHdx3hrx/dw6O7jlDwGOg2Ayox1tAaYHfJ/B7g4knqXwd8vWS+Q1IvxQfb3xoRXx1vJUmbgc0AZ555ZlkBm9UCj29ks2VWO4slvQ3oAf6opHh9RPQAbwU+Luns8daNiC0R0RMRPcuXL5+FaM2qy+Mb2WypRCLYC6wrmV+blJ1E0mXAh4A3RsTIaHlE7E3enwIeAC6qQExmNc/jG9lsqUQieBjYKOksSW3AVcBJV/9Iugj4C4pJ4FBJ+WJJ7cn0MuBVQGkns1nD8vhGNlvKTgQRkQOuB+4FngDuiojHJd0i6Y1JtT8CuoAvj7lM9FygV9L3gW9T7CNwIjDD4xvZ7KnIg2ki4h7gnjFlN5VMXzbBev8CXFCJGMzqzej4Rr/63t8k37WS69+12U9FsxnhO4vN5jCPb2SzwYnAzKzBORGYmTU4JwIzswbnRGBm1uAqctWQmdWPQiHYtnuAZ/qPs2HpfF+p1ACcCMzsBI9v1JjcNGRmJ3h8o8bkRGBmJ3h8o8bkRGBmJ3h8o8bkRGBmJ3h8o8bkRGBmJ4yOb9S1/avMe/ofec+lG91R3ACcCMzsJB7fqPE4EZiZNTgnAjOzBucbysys6nw3c3U5EZhZVflu5uqrSNOQpMslPSlpp6QbxlneLunOZPlDkjaULLsxKX9S0usqEY+Z1Q7fzVx9ZScCSc3AJ4DXA+cBV0s6b0y164AjEXEOcBvwB8m651F82P35wOXA/0k+z8wahO9mrj5FRHkfIL0SuDkiXpfM3wgQEb9fUufepM53JLUAB4DlwA2ldUvrTbbNJevPjdd88PZpx7rt+9sA2HThpmmvm8kVSOfy015vx/bHANh43kunve7p8ja9zVra5tF0jr0DKUoPRRKs6Z7Hgg63Xpea19pMa/Pp/36/610//UhE9Iwtr8S3vAbYXTK/B7h4ojoRkZM0CCxNyr87Zt01421E0mZgM0DXqrNPK9DTSQDlms3/vN6mt1mL2+xqb2ZeazOpbJ6IYhKY19pMV/vMNw7Ue5KdqppJtxGxBdgC0NPTE3f+6itndft7B1I82z88q9s0axTVumro3X95IwA3vX/rjG+rEtt80fL5rFzYcdrbvutd45dXIhHsBdaVzK9NysarsydpGloE9E9xXTOrc01N4mXrF/Oy9YurHUpDqsRVQw8DGyWdJamNYufv2FS3FbgmmX4T8K0odk5sBa5Krio6C9gI/GsFYjIzsykq+4wgafO/HrgXaAZuj4jHJd0C9EbEVuBTwOck7QQOU0wWJPXuArYDOeDdETH9HlkzMzttFekjiIh7gHvGlN1UMp0GfmmCdT8KfLQScZiZTVWhEGSWnkO+ayWP7jrS0Hcz10xnsZlZpYzezXzsvCuhuYU//daOhr6b2YPOmVnDGb2bmZY2UFPD383sRGBmDcd3M5/MicDMGo6fzXwyJwIzazijz2Zub2lC0PDPZnZnsZk1nNFnM/sZCEVOBGbWkKpxN/NcvWTVicDMbBbM5UtW3UdgZjYL5vIlq04EZmazYC5fsupEYGY2C+byJatOBGZms2AuX7LqzmIzs1kwly9ZdSIwM5slc/UBPG4aMjNrcE4EZmYNzonAzKzBlZUIJC2RdJ+kHcn7Cxq+JG2S9B1Jj0v6gaS3lCz7jKSnJW1LXpvKicfMzKav3DOCG4D7I2IjcH8yP9Yw8PaIOB+4HPi4pNLrpX47IjYlr21lxmNmZtNUbiK4Argjmb4DuHJshYj4UUTsSKb3AYeA5WVu18zMKqTcRLAyIvYn0weAlZNVlvRyoA34cUnxR5Mmo9sktU+y7mZJvZJ6+/r6ygzbzMxGnTIRSPqmpMfGeV1RWi8iAohJPmcV8DngHRExOuDGjcBLgJ8ClgAfmGj9iNgSET0R0bN8uU8ozMwq5ZQ3lEXEZRMtk3RQ0qqI2J8c6A9NUG8h8HfAhyLiuyWfPXo2MSLp08BvTSt6MzMrW7lNQ1uBa5Lpa4Cvja0gqQ34CvDZiLh7zLJVybso9i88VmY8ZmY2TeUmgluB10jaAVyWzCOpR9InkzpvBn4OuHacy0S/IOmHwA+BZcBHyozHzMymqayxhiKiH3j1OOW9wDuT6c8Dn59g/UvL2b6ZmZXPdxabmTU4JwIzswbnRGBm1uCcCMzMGpwTgZlZg3MiMDNrcE4EZmYNzonAzKzBORGYmTU4J4Ip6mproa1F1Q7DzKzinAimaFFnKxeu7WZN9zyanA/MrI44EUxDS3MTZy7t5MJ13Szraqt2OGZmFeFEcBo6WpvZuHIBL12zkAUdZY3bZ2ZWdU4EZVjQ0cpL1yziJ1Z20d7qr9LMapN/zlbA0q52Fne2cWAozd6BFLn8hE/sNDObc/wztkKamsTq7nlsWtfNGYs6kDuUzaxGOBFUWGtzE2ctm8+Fa7tZPL+12uGYmZ1SWYlA0hJJ90nakbwvnqBevuQxlVtLys+S9JCknZLuTJ5vXBfmtTXzkjMWct6qhcxvb652OGZmEyr3jOAG4P6I2Ajcn8yPJxURm5LXG0vK/wC4LSLOAY4A15UZz5yzqLOVC9Ys4uzl831DmpnNSeUmgiuAO5LpO4Arp7qiJAGXAnefzvq1RBIrFnawad1i1i6eR7PvSDOzOaTcRLAyIvYn0weAlRPU65DUK+m7kkYP9kuBgYjIJfN7gDUTbUjS5uQzevv6+soMuzqam8S6JZ1cuG4Ryxe0VzscMzNgCpePSvomcMY4iz5UOhMRIWmi6ybXR8ReSS8CviXph8DgdAKNiC3AFoCenp6avj6zvaWZc1Z0ccaiDnb1H2colTv1SmZmM+SUiSAiLptomaSDklZFxH5Jq4BDE3zG3uT9KUkPABcBfwV0S2pJzgrWAntPYx9qVld7C+evXsTh4xl29R8nnS1UOyQza0DlNg1tBa5Jpq8Bvja2gqTFktqT6WXAq4DtERHAt4E3TbZ+I1gyv41N67rZsKyTlmb3H5jZ7Co3EdwKvEbSDuCyZB5JPZI+mdQ5F+iV9H2KB/5bI2J7suwDwG9I2kmxz+BTZcZTsySxatE8LlrXzSrfkGZms0jFH+a1paenJ3p7e6sdxoxKZ/Ps6h/m8PFMtUMxszniRcvns3Jhx2mvL+mRiOgZW+47i+eojtZmXnzGAs5bvZCudg8JZWYzx4lgjls0r5UL1i7inBVdtLX4z2VmleefmjVi+YJ2ls5vY99gin0DafKF2mvSM7O5yT8xa0hTk1i7uJNN67pZsbDdHcpmVhFOBDWoraWJs5d3ccGaRSya5xFOzaw8TgQ1bH57C+etXshLzljAvDaPcGpmp8d9BHVg8fw2ujtbOTg0wp4jw2T9hDQzmwYngjohiTMWdbCsq419A2n2D6Zwf7KZTYUTQZ1paW7izKWdrOru4OBQmoNDaTI5ZwQzm5gTQZ1qbW5i7eJOVi+ax3PHRtg/mGY4k692WGY2BzkR1LmmpuJDcVYs7GBgOMO+gTSDqWy1wzKzOcSJoIF0d7bR3dnGcCbHvoE0/cdG3I9gZk4EjaizrYVzVnRx5pLOE/0IvtLIrHE5ETSwtpYm1i3pZHV3sR9h30DKD8cxa0BOBEZzk1i5sIMVC9o5Mpxl/2DKj880ayBOBHaCJJbMb2PJ/DaOjeTYP5Ci/3iGGnxkhZlNQ1lDTEhaIuk+STuS98Xj1Pl5SdtKXmlJVybLPiPp6ZJlm8qJxyqnq72FjSsXcNGZ3azu7vAjNM3qWLljDd0A3B8RG4H7k/mTRMS3I2JTRGwCLgWGgW+UVPnt0eURsa3MeKzC2luaWb90Pi87czEblnXS3urhqczqTbn/q68A7kim7wCuPEX9NwFfj4jhMrdrs6y56flnKv/Eyi4WdLhV0axelJsIVkbE/mT6ALDyFPWvAr44puyjkn4g6TZJ7ROtKGmzpF5JvX19fWWEbOWQxNKudl66ZhHnr1nI0q42PxfBrMad8medpG8CZ4yz6EOlMxERkibsVpS0CrgAuLek+EaKCaQN2AJ8ALhlvPUjYktSh56eHndfzgELO1pZ2NFKOptn/2DxBjXfj2BWe06ZCCLisomWSTooaVVE7E8O9Icm+ag3A1+JiBPjG5ScTYxI+jTwW1OM2+aQjtZmzlo2nw1LOzmeyTMwnGEwleVoOucrjsxqQLkNvVuBa4Bbk/evTVL3aopnACeUJBFR7F94rMx4rIok0dXeQld7C2sXQ74QDKayDKayDAxnfLOa2RxVbiK4FbhL0nXALoq/+pHUA7wrIt6ZzG8A1gH/MGb9L0haDgjYBryrzHhsDmluev6+BJhPOps/kRgGU1lybkYymxPKSgQR0Q+8epzyXuCdJfPPAGvGqXdpOdu32tLR2kxHazMrF3YQERwdyTE4XEwKx0bcjGRWLb4G0KpC0onO5nVALl94vhkplWXEzUhms8aJwOaEluYmlna1s7SreAVxKpNPkkKGoVSOvMfLNpsxTgQ2J81ra2ZeWzNnLCo2Iw2lcwylsgwMZzmecTOSWSU5EdicJ4lF81pZNK+VdUsgmzQjDQxnOT6SI53N+wE7ZmVwIrCa09rcxLKudpYlzUgRwUiuQCqTJ5VNXpk86WzeN7iZTYETgdU8SSeuSBo7/G02XyCVzZPO5ElnCycSRTqbd/OSWcKJwOpaa3MTrc1NLOxoPak8Ik5KDKNnEKls3vc3WMNxIrCGJOlEh/RYmVzhxFlDKpMnnSu+j+QKPouwuuREYDZGW0sTbS1NLJp38llEoRCkc3lGsgWyhQK5fJDLx4npbL5ANl8gVwifVVhNcSIwm6KmJtHZ1kJn26nrRgTZfJArFMgmSWI0WRQTRaFkeXHaZxtWLU4EZjNAEm0tom0aj/zIJUniRNJIkkguXyBfCAoBhYjiq1CcjqQsH0FEUiepazZVTgRmc0RLcxMtzcUxmSqhmBBGk8RoEkkSSOH5BFIYu7zw/DxAUDzDGZ0GSs5enj+TGbssSpeVlJ28/guXzbbJzsQmiygmWXHy9YrtWIsAAASVSURBVCaNZtL1ZuohUE4EZnWqqUk0If8nt1Pyk8jNzBqcE4GZWYNzIjAza3BOBGZmDa6sRCDplyQ9LqmQPJ5yonqXS3pS0k5JN5SUnyXpoaT8TklTuELbzMwqqdwzgseAXwQenKiCpGbgE8DrgfOAqyWdlyz+A+C2iDgHOAJcV2Y8ZmY2TWUlgoh4IiKePEW1lwM7I+KpiMgAXwKukCTgUuDupN4dwJXlxGNmZtM3G30Ea4DdJfN7krKlwEBE5MaUj0vSZkm9knr7+vpmLFgzs0ZzyntNJH0TOGOcRR+KiK9VPqTxRcQWYAtAT0+Pb6A3M6uQUyaCiLiszG3sBdaVzK9NyvqBbkktyVnBaPkpPfLII89J2nWa8SwDnjvNdWuJ97O+eD/rS7X2c/14hbNx9/nDwEZJZ1E80F8FvDUiQtK3gTdR7De4BpjSGUZELD/dYCT1RsSEVzjVC+9nffF+1pe5tp/lXj76C5L2AK8E/k7SvUn5akn3ACS/9q8H7gWeAO6KiMeTj/gA8BuSdlLsM/hUOfGYmdn0lXVGEBFfAb4yTvk+4A0l8/cA94xT7ymKVxWZmVmVNOKdxVuqHcAs8X7WF+9nfZlT+6nJxtQ2M7P614hnBGZmVsKJwMyswTVUIpho8Lt6ImmdpG9L2p4MCPjeasc0UyQ1S/qepL+tdiwzSVK3pLsl/X9JT0h6ZbVjqjRJ70v+vT4m6YuSOqodU6VIul3SIUmPlZQtkXSfpB3J++JqxtgwieAUg9/VkxzwmxFxHvAK4N11up8A76V4SXK9+xPg7yPiJcCF1Nk+S1oDvAfoiYiXAs0U7zeqF58BLh9TdgNwf0RsBO5P5qumYRIBEwx+V+WYKi4i9kfEo8n0UYoHjQnHcKpVktYC/wH4ZLVjmUmSFgE/R3KPTURkImKgulHNiBZgnqQWoBPYV+V4KiYiHgQOjym+guJAmzAHBtxspEQw0eB3dUvSBuAi4KHqRjIjPg68HyhUO5AZdhbQB3w6aQb7pKT51Q6qkiJiL/A/gWeB/cBgRHyjulHNuJURsT+ZPgCsrGYwjZQIGoqkLuCvgF+PiKFqx1NJkv4jcCgiHql2LLOgBXgZ8OcRcRFwnCo3I1Ra0j5+BcWktxqYL+lt1Y1q9kTxGv6qXsffSIlgosHv6o6kVopJ4AsR8dfVjmcGvAp4o6RnKDbxXSrp89UNacbsAfZExOhZ3d0UE0M9uQx4OiL6IiIL/DXw01WOaaYdlLQKIHk/VM1gGikRnBj8Lnkk5lXA1irHVHHJA38+BTwREf+r2vHMhIi4MSLWRsQGin/Hb0VEXf6CjIgDwG5JL06KXg1sr2JIM+FZ4BWSOpN/v6+mzjrEx7GV4kCbMI0BN2fKbIw+OidERE7S6OB3zcDtJYPf1ZNXAb8M/FDStqTsg8l4T1abfg34QvID5ingHVWOp6Ii4iFJdwOPUrzq7XvMsSEYyiHpi8AlwLJkkM7fBW4F7pJ0HbALeHP1IvQQE2ZmDa+RmobMzGwcTgRmZg3OicDMrME5EZiZNTgnAjOzBudEYGbW4JwIzMwa3L8BpvJKhJkZ07cAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from statsmodels.graphics.tsaplots import plot_acf, plot_pacf\n",
    "import pandas as pd\n",
    "\n",
    "pdata=pd.read_csv(\"http://image.cador.cn/data/airmiles.csv\")\n",
    "plot_acf(pdata.miles, lags=10, title=\"airmiles autocorrelation\")"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
